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Abstract 



0^ , High spin states of ^^S and ^^Ni are investigated by means of the cranking Hartree- 

I Fock method with the Gogny interaction without imposing a restriction on the axial 

' reflection symmetry. It was found that a non-axial octupole deformation of the Y^i type 

■ becomes important in the yrast states of ^^5. A similar effect is predicted for the nucleus 

g : ^^Ni. 

PACS number(s): 21.60-n; 21.60.Jz; 27.30.+t; 27.40. +z 

§ ; 1 Introduction 

One of the most important concepts in the many-body theory of finite Fermi systems is the 
^ ' mean field approach. In fact, many phenomena observed in nuclei can be explained by means of 
' a spontaneous symmetry breaking mechanism which leads to a mean field solution that does not 
obey symmetries of the original many-body Hamiltonian Q . A nucleus with a partially filled 
highest shell spontaneously deforms in such a way that a lower energy minimum is achieved, and 
thus a non-spherical equilibrium shape is attained. Coriolis forces may cause a similar effect 
for a nucleus which has a spherical shape in a ground state. In particular, a strong coupling 
of normal and intruder states near the Fermi surface at large rotational frequencies can lead 
to a superdeformed or to an octupole deformed shape for certain combinations of protons and 
neutrons |^ . The most practical method for the analysis of nuclear shapes is a phenomenological 
macroscopic-microscopic method (MMM) which combines the liquid drop model describing the 
bulk properties of nuclear matter and the Strutinsky Shell Correction method providing the 
description of quantum shell effects of phenomenological potentials . The more fundamental 
approach is based on self-consistent Hartree-Fock (HF) calculations once a particular choice 
of nucleon-nucleon interaction has been made. A commonly used inter-nucleon force is that 
of Skyrme. While the HF+Skyrme approach describes major features of nuclei quite well, it 
does not completely take into account the pairing correlations. In addition, various sets of 
parameters for Skyrme forces may not provide a definite answer in some cases. The Hartree- 
Fock-Bogoluibov (HFB) method with the Gogny forces resolves these problems quite effectively. 
Furthermore, it has the same prediction power as the MMM 0]. 
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In recent years a lot of experimental and theoretical efforts has been devoted to the analysis 
of superdeformed (SD) bands in different mass region ^, ^, On the other hand, the 

study of octupole degrees of freedom is also a topical subject in the nuclear structure physics 
1^, It turns out that octupole deformations are significant for superdeformed nuclei as 
well as for normal-deformed (ND) nuclei. Most of these studies were restricted to the axial 
octupole deformation proportional to the Y30 term in the octupole family. The reason for this 
is primarily computational complications arising as a result of the extra degree of freedom 
introduced by non-axial octupole deformations. However, during the last few years some re- 
markable results have been reported. It was found that Y31 and 132 components resulted in the 
lowest energy octupole vibrations for oblate and prolate superdeformed shapes and therefore 
may have consequences for octupole instability |Tl], |^. Calculations using the MMM with 



the Woods-Saxon potential predicted the importance of the banana-type Y31 deformation for 
highly deformed nuclei [|l3l, while the Y32 deformation has been found to be important in the 
^^"^Ra nucleus The self-consistent HF+Skyrme calculations |15] suggest the softness of 

the oblate states in A ~ 80 nuclei against the 133 deformation in the ground states. The self- 



consistent cranking HF-|-Skyrme approach |T6| predicts that the 131 deformation is important 
for a correct description of the yrast band of ^"^S. Moreover, the study of octupole deformations 
sheds light on the interplay between a classical chaos and a quantum spectrum of finite Fermi 
systems |17, 18, 19 1. Guided by thorough investigation of non-axial octupole deformations in 
a harmonic oscillator model [^, we aim to analyze how symmetries break at high spins in 
the cranking HFB approach with the Gogny interaction. In accordance with the octupole in- 
stability suggested at particle number =16 and 28 |20|, we choose two N = Z nuclei, ^"^S 
and ^^Ni. In Sec.^ we review the main features of our model. The discussion of main results 
is presented in Sec.| with the conclusion following in Sec. 4. In the Appendix we introduce a 
simple two-level model to understand qualitatively the behavior of the Inglis-Beliaev moment 
of inertia due to strong octupole coupling between two single-particle (s.p.) states. 



2 The Model 

As mentioned above, the present self-consistent cranking HFB calculation has been performed 
with the effective Gogny DIS interaction |21|, |22|, p3|, E4l. This interaction provides a good de- 



scription of many nuclear properties over the nuclide chart, e.g., ground state energies, odd-even 



energy differences, electron scattering data ^2[, fission barriers Also, a good description is 
obtained for bulk properties of rotating nuclei in actinide region p5|, mercury region ^6] and 
fp-shell region . Our numerical code has been used for the analysis of microscopic dynamics 



in light rotating nuclei [^, |3^ . In these calculations, however, the signature symmetry was 
conserved (see discussion about different symmetries in rotating nuclei in ^). The present 
calculations are performed without assuming d priori the axial and signature symmetries. In 
addition, our Hamiltonian includes the Coulomb interaction and the center of mass correction 
up to the exchange terms exactly. 

The 2;-axis is taken as a rotational axis in our code. To save the CPU time in numerical 
calculations, we impose the Pe"*'^"^^ (z-simplex) and the Pe~'"^'^^f (Sy) symmetries where 

P is the parity operator, e"*'^"^' is the rotation operator around the i-axis (i=y,z) on angle 
TT , and f is the time-reversal operator. Due to the z-simplex and Sy symmetries, the mass 
asymmetry of nucleus is allowed only along the x-axis. Thus, we solve numerically the following 
cranking HFB equations: 



S{(piuj)\H - XpZ - XnN - LOJ, 
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+^l,{(t){Lo)\x\<P{u))x\(t){uJ)) = 0. (la) 
(0(c.)|Z|0(c.)) = Z, (0(c.)|iV|0H) = N, (lb) 
(0(c.)|J,|0(c.)) = I, (Ic) 

(0(c.)|x|0(c.)) = 0, (Id) 

where the Lagrange multiphers Xp and A„ are the chemical potentials of proton and neutron 
fields, respectively (the operators Z and are protons and neutrons number operators); the 
Lagrange multiplier u is the angular frequency of a collective rotation around the 2;-axis and 
is the 2;-component of the angular momentum operator J. To keep the center of mass motion 



fixed, we also impose the quadrupole constraint operator ij,x{(j){oj)\x\(j){uj))x ||32| to the Routhian 
R = H ~ XpZ — A„A^ — ujJz in the x-axis direction. 
The quasi-particle (q.p.) operators 



= ^kici + VkiCk , (2) 

k 

where the state \k) is a single-particle basis state (see below), obey the equation of motion 

[R,a^i\ = ei{uj)a\ (3) 

which defines the quasi-particle energies eiiuj) and quasi-particle amplitudes Uki and Vki as 
functions of the rotational frequency uj. 

We define the y-axis as the quantization axis of deformation for convenience. Consequently, 
the quadrupole, [32 and 7, and octupole , /Ssm, deformation parameters are defined as 

= T ARl 

F((3y2-r2)) 



5 ARi 

5 Y IARq 

[3^ {{x^-z^)) 
" V 5 ARl ' 

= T ARl 

V 9 ARl 

n _ 47r (r3(y3i(g,y)-y3„i(g,y))) 



'2l7r (x(5?/2 - r^)) 
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ARl 



n ^ 47r {r\Y,,{d,^)-Y,_,{d,^))) 



/357r - 3^2)) 

where (r, 6, ip) are spherical coordinates related to the Cartesian coordinates in the rotating 
frame {x,y,z) as 

{x,y,z) = {r sin 6 cos !f , r cos 6 , r sin 6 simp) . (10) 

All deformations proportional to the odd power of the y coordinate are forbidden in our code 
due to the Sy symmetry. Since /?3o and /332 are proportional to the odd y"' terms, we use /?3i 
and to represent the degree of the non-axial octupole deformation, when the y-axis is the 
largest axis of a prolate nucleus. Notice that the octupole deformation with a mass asymmetry 
like a pear shape is also represented by the combination of /?3i and when the x-axis is the 
largest axis of a prolate nucleus Q. 

The s.p. wave functions have been expanded in a three-dimensional harmonic oscillator 
basis up to the principal quantum number = 8 for ^^5* and up to = 10 for ^^Ni. The 
basis has been symmetrized with respect to the z-simplex operation, and eigenfunctions are 



eigenstates of the Sy operator (the Goodman basis ||33|). Since the ground state shapes of the 
chosen nuclei are a normal-deformed one for ^"^S and a spherical one for ^^Ni, we use a spherical 
Cartesian basis with the same range parameters of the Hermite polynomials for all axes. The 
range parameters have been optimized to reproduce the largest binding energy of each ground 
state. 

To understand the relation between s.p. degrees of freedom and collective effects due to the 
rotation, we calculate three moments of inertia, the static moment of inertia J'^^\lo) = I /lu, the 
dynamical moment of inertia J'^'^\uj) = dl /du, and the Inglis-Beliaev (IB) moment of inertia 



\Jij{' 



j,m=2j: ;7:"\ (12) 

i>j + ej(i^) 

which is the leading order term of the nuclear moment of inertia. Here Jij{uj) is a matrix element 
of the angular momentum operator 

Jij{uj) = {(l){uj)\[ajai, Jj|0(w)). (13) 
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It is well known that the IB moment of inertia is too small to reproduce the absolute value 



of the moment of inertia. Comparing JTibIw) and J''^^\uj) or J'^'^\lli), we can find the validity of 
^From these /?3i and P33 values, for an example, we can obtain the value of for the x-quantization axis 

as 

/33o = ^/33r- 4^/333. (11) 
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using the J'mi^^) for the analysis of structural changes in the nuclear mean field. In addition, 
we analyze different contributions to the IB moment of inertia, Eq.(|T2|), 

i7ib(w) = ^ Jt,s{'^) = Jp,+li^) + Jp-l{^) + i7n,+l(a;) + Jn-l{uj), (14) 

T,S 

where (r, s) denote the isospin t = p,n and the z-simplex quantum number s = ±1 which char- 
acterize different subspaces. The comparison of different components of Jisi^^) could provide 
the information in which subspace the q.p. degrees of freedom mainly affect the bulk properties 
of the nuclei. 

The cranking HFB equations ([Ta|) are solved in an iterative way. As the convergence con- 
dition for each HFB state, we impose the condition 

< .l[KeV], (15) 
where n is the number of iterations in the course of solving the cranking HFB equations. 



E(n) , X (n— 1), X 

\4 (^)-ei (^) 



3 Discussion of results 
3.1 

In Figs. 1-3 results for the yrast line states, angular momenta and deformations as functions 
of the rotational frequency are presented. The binding energy of the ground state in ^'^S is 
well reproduced in our calculations (Fig.|I]). We found that the contribution from the pairing 
interaction terms to the total binding energy is almost negligible. The energy gaps between 
q.p. energies near the Fermi energies are large, about 4 MeV, both in the proton and neutron 
fields. Therefore, the calculations for the ND and SD bands have been done using the cranking 
HF approximation with the same parameter set as for the HFB calculations. While there are 



differences in details between our calculations and those with Skyrme forces ||I6[, we obtained 
similar main results for the ND and SD bands. The SD band becomes the yrast one for / > 12h. 
The excitation energy of the SD minimum relative to the ground state is about 10 MeV (Fig. 
|ip. Except for the binding energy, these results are also consistent with those of SLy4-HF 
calculations ||3^, in which /5-deformations are (3 = .16 and (3 = .7 in the ground state and 



SD minimum state, respectively. In comparison with the results of the Skyrme HI calculations 
T6| , however, finite values of the non-axial octupole deformation Y31 in the ND yrast band are 



obtained at larger rotational frequency u > 1.5 [MeV/ h]. As is seen in Fig.^(a), the value 
of the non-axial octupole deformation |/53i| increases suddenly at u; ~ 1.5 [MeV/ h], where 
the pseudo level-crossing occurs between s.p. orbits in the subspaces with s = ±1 (Fig.|^(a)), 
both in the proton and neutron fields. These orbits can be associated with the principal 
quantum numbers N = 2 and N = 3 at uj = 0[MeV/h]. As a result, the Y^i deformation 
becomes favorable at uj ~ 2.2 [MeV/ h], for ^^5* with a finite value /Ssi = .13. The density 
distribution projected on the plane perpendicular to the rotation axis is shown on Fig.|^(a) for 
the octupole band at a; ~ 2.2 [MeV/^] where the jS^i deformation has a maximal value. High 
resolution 7— spectroscopy could test the validity of different predictions based on different 
nuclear interactions. 

From the analysis of Figs.^(a), |^(a) it follows that there is a correlation between the behavior 
of the dynamical moment of inertia J'^'^^co) and the occurrence of the octupole deformation. 
Similar to the behavior of the |/?3i|, the dynamical moment of inertia J'^'^\uj) begins to increase 
at a; ~ 1.5 [MeV//i], and both quantities show local maxima around a; ~ 2 [MeV/^]. Due 



5 



to the strong octupole interaction the quasi-crossing between s.p. levels is smooth (Fig.|^(a)). 
The static moment of inertia J''^^\oj) is less sensitive to the quasi-crossing of s.p. levels, and it 
changes smoothly with the increase of the rotational frequency uj. However, the IB moment of 
inertia may reflect the structural changes in a mean field like the dynamical moment of inertia 
J'^'^\uj). In fact, it also begins to increase at uj ~ 1.5 [MeV//i] and shows a local maximum 
at uo ~ 2.2 [MeV/^]. The correlation between the IB moment of inertia and the magnitude of 
the octupole deformation |/33i| may be understood within a simple two-level model presented 
in the Appendix 0. According to this model, the two rotation-aligned single-particle states 
with a strong octupole coupling near the Fermi-surface may give a large contribution to the 
value of the IB moment of inertia. From Fig.|^(a) follows that the magnitudes of Jp^+i{uj) 
and J7n,+i(cj) (see Eq.(|14[)) are smaller than those from the s = — 1 subspaces which play the 
dominant roles in the collective rotation at low spins, in the region uj < 1.5[MeV/^]. However, 
at the UJ > 1.5[MeV/^], and J7n,+i(t^) begin to increase rapidly, and they show local 

maxima at ci; ~ 2.2[MeV/^] where also shows the maximum. Since J^^_i(cj) and J'n-i{cd) 
are almost unchanged in this region, the occurrence of the Y31 deformation is mainly due to 
the s.p. quasi-crossing in the s = +1 subspaces. To justify this statement we calculate the 
quantity cos'^ipi^) (see Appendix). Using the identity of the intruder orbits /i which have odd 
parity at = 0[MeV//;.], we evaluate the cos^ ip{uj) approximately as: 

cos^VC'^)- \{n^,ny,n^;a\fiiu;))\'^, (16) 

nx +ny +nz =odd 

where a s.p. vector jn^,, n^, n^; a) is a component of the Goodman basis in the s = +1 subspace, 
and 71x1 ny and are the number of quanta on the x— , y— and z— axes, respectively. As 
is shown in Fig.^(a), the evaluated value of cos'^ ip{uo) = 1/2 or equivalently {ipiuj)] = 7r/4 at 
u ~ 2.2 [MeV/^]. According to the analysis within the two-level model (see Appendix), at this 
value of \ipiuj)\ the mixing of s.p. states with opposite parities is the strongest one. 

Similar to the results of []TB|, our calculations show rather shallow minima for the jSsi 



deformation. In Fig.|TO|(a) we present results for the fluctuations of octupole and quadrupole 
deformations 



(see also Fig.|ll](a)). The contribution of the ground state octupole correlations may improve 
the mean field results. This contribution could be estimated, for example, in the random phase 



approximation like according to the prescription [^, however, this is beyond the scope of 
the present paper. 

As is seen in Fig.|^(a), both moments of inertia, J'^'^\uj) and J/ib(w), begin to increase again 
at a; ~ 2.5 [MeV/^] , which is apparently due to the quasi-crossings in the s = — 1 subspaces 
at this point. There are also level-crossings between s.p. orbits with different simplex numbers 
(Fig.^a)), which suggest the instability of the 2;-simplex symmetrized state at high spins. We 
may speculate that the breaking of the z-simplex symmetry due to the tilted or the chiral 
rotation |jlO| may lead to a better description of high spin states in ^^5* at uj > 2.5 [MeV/Zi]. 



3.2 ^^Ni 

Our calculations reproduce the binding energies and the charge radii of ground states in se-eojy^ 
quite well (Fig. |T2|). In this chain of nuclei we concentrate our attention on ^^Ni. 

Though the shape of the ground state in ^^Ni is spherical, we also found two shape isomers 
with the ND and SD configuration at u; = [MeV/^], which are ~ 5 MeV and ~ 18 MeV 
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above the ground state, respectively (see Fig.^). The yrast non- rotating states observed in 
^^Ni are related to vibrational excitations. In fact, their properties can be studied within 



the self-consistent cranking approach with the rotation around a symmetry axis and we 
will discuss this subject in a forthcoming paper. The calculations for the ground state, ND 
and SD bands have been done in a similar manner to those of ^^S* in the HF approximation, 
since the contribution of the pairing interaction terms to the total binding energies are almost 
negligible. On Fig.|^ the evolution of the angular momentum with the increase of the rotational 
frequency is shown for SD and ND bands. At low spins the SD minimum is formed by the Ap-4h 
configuration with respect to the ground state in ^^Ni, 7r[(3)~^(4)^]z/[(3)~^(4)^], which has a 
prolate shape with the quadrupole deformation /? = .6 (Fig.|^). This is also Ah configuration 
with respect to the SD state in ^°Zn due to the N = 30 SD gap in the Zn isotopes [H5 



which is also seen in Fig.|Tj. The MMM calculations ^ also predict the SD shape at large 
rotational frequencies. We restricted, however, our analysis of the SD band in the region of 
uj = 0.0 — 0.4[MeV/^], since at larger rotational frequencies the contribution of higher shells 
(A^ > 10) becomes important. 

Using the decomposition of the rational harmonic oscillator (RHO) into the isotropic ones. 



the relation between multi-clusters and mean field results has been discussed in |^8[. For SD 
shapes in the RHO the one symmetric and one asymmetric combinations of spherical oscillators 
are expected (see also |12|). Within the RHO the SD states in ^^Ni correspond to the asym- 



metric combination of two spherical oscillators with magic numbers 40 and 16. The density 
distribution of the SD state (see Fig.|15D does not show such a multi-cluster structure at = 
[MeV/^]. The possible reason is that our configuration space is too small. However, the SD 
minimum could be related to the resonance state in the ^^Si+^^Si collision at high excitation 



energy -E^m = 65 — 70MeV and high angular momenta J = 34 — 42/i reported in [S^, A 
thorough study of ^^Ni as well as ^'^S at high spins may help to understand the link between 
the tendency for nuclei to create strongly deformed shapes and the tendency to develop the 
cluster structure. 

A prolate deformation f5 = .35 is found for the ND minimum (Fig. ^ at a; = 0.0 [MeV//i]. 
With the increase of the rotational frequency the non-axial octupole deformation Y31 arises in 
the ND band. As is seen in Fig.^(b), the value of the non-axial octupole deformation increases 
rapidly around uj ~ .9[MeV//;,]. It corresponds to / ~ 10 [/?.], where the pseudo level-crossing 
occurs between s.p. orbits which have the principal quantum numbers = 3 and A^ = 4, 
respectively, at u; = 0[MeV/^] (Fig.|](b)). The maximal value of the F31 deformation, jS^i = .09, 
is approached at c<j ~ 1.4[MeV/^] (Fig.^b)). As is seen in Fig. ||(b), J^,+i(a;) and J7n,+i(i^) 
increase rapidly around u ~ .9[MeV/^] and approach the maximal value at ~ 1.3 [Me V/^]. 
Using the two-level model and similar arguments as for the case of ^^S* we may conclude: since 
J'p^_i{lj) and J7„,_i(a;) are decreasing in this region of the rotational frequency, it is most likely 
that the octupole deformation is caused by the quasi-crossing of s.p. levels, which have different 
parities ai u = 0[MeV/h], in the s = -|-1 subspaces. Both cos^ ip{uj) of the intruder orbits with 
the positive parity at a; = 0[MeV//;.] in the proton and neutron fields. 



cos 



rix+ny+n 



have maximal values at ~ 1.3[MeV//;,] and u ~ 1.35[MeV//i], where the components of the 
IB moment of inertia, J'p,s=+ii^) and J7n,s=+i(w), show maxima, respectively (Fig. || (b)). 

The dynamical and IB moments of inertia also increase at ~ .9 [MeV/^] and show maxima 
around to ~ 1.2 [MeV//i] (Fig.|^(b)). However, the behavior of the J'^'^\uj) and J^ib{uj) moment 
of inertia is less correlated at larger values of the rotational frequency. The pseudo-crossing is 
much sharper in ^^Ni in comparison with the one in ^"^S (Fig.^), i.e. the octupole interaction 
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is expected to be weaker. Both cos^^/'(w) decrease rapidly at > 1.55[MeV/^] (see Fig.(^b)). 
Due to the termination of the parity mixing of the s.p. levels, the matrix element Jij{uj) between 
mixed states, Eg. ([T3|) , becomes rapidly small. As a result, the IB moment of inertia, Jisi^^), 
also decreases. On the other hand, the dynamical moment of inertia, J'^'^\u;), reflects changes 
of the self-consistent mean field and, in particular, the modification of the two-body interaction 
due to the rotation. Due to the two-body interaction, the high intruder unoccupied state 
[s = +) contributes to the sudden increase of the dynamical moment of inertia observed at 
high rotational frequencies u > 1.5[MeV/h]. 

The fluctuations of octupole deformations AQsi and AQ^s are always larger than 1 in all 
region of the ND band (Fig. |lO|(b)), in spite of the presence of the octupole minimum (Fig. 
|TT](b)). The octupole minimum in ^^Ni is much more shallow in comparison with the one of 
^"^S and it could explain large fluctuations of the octupole deformations in the considered case. 

4 Summary 

It was expected from various calculations based on the MMM that octupole deformations 
would arise in rotating nuclei 0. Using the cranking HF(B) approach with the effective Gogny 
interaction, we found that the non-axial octupole deformation associated with the Y31 term 
in the octupole family becomes important in the yrast band of ^"^S at the angular momenta 
I > 5[h]. The primary mechanism behind the occurrence of the octupole deformation is related 
to the strong mixing via octupole interaction of s.p. orbits with a positive simplex quantum 
number. Similar phenomena have been observed in the nucleus ^^Ni, where we predict the 
octupole softness in the excited ND band at high spins / > 10[h]. 

Finally, an exploration of the octupole phenomenon certainly could deepen our understand- 
ing of different aspects of a spontaneous symmetry breaking mechanism in finite Fermi systems 
like nuclei. In particular, the breaking of the intrinsic reflection symmetry could be related to 
unexpected strong electric dipole and octupole transitons in rotational bands and to a formation 
of multi-cluster structures in strongly deformed nuclei. Measurements using new generations 
of modern detectors can test the predictions made within the HF(B) approach and lead to new 
insights regarding effective nucleon-nucleon interactions and their properties. 
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A A two- level model 

Let us consider the situation such that the high-lying intruder s.p. orbit n is coming down to 
the highest-lying occupied orbit u due to the Coriolis term —ujJz] they have the same simplex 
quantum number. For the sake of simplicity, we restrict our discussion within the 2x2 subspace. 
We assume that i) our states are eigenstates of the parity operator P at = 0; ii) there are no 
changes in the subspace expanded by the simplex partners of /i and z/ at 7^ 0. 
At Q < uo < uOc the eigenstates can be written as 

|i/(ci;)) = COS'?/'(w)|z/(tj = 0)) + sin'?/'(a;)|/i(tj = 0)), 
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\fi{uj)) = cos^/'(w)|/i(w = 0)) - sin?/'(w)|z/(w = 0)), (A.l) 
where is a monotonic function of u 0, which satisfies the following conditions; 

0<|^/'H|<|, (A.2) 

V'(c^ = 0) = 0, \^{u,)\ = ^. 

Here, Uc is the largest rotational frequency. Using the unitary transformation in Eq. (|A.l|) , at 
u ^ the density matrix has the following form in this subspace 



_ I cosijj{uj) sinipiuj) \ / 1 \ 
\^ —smip{uj) cosiP{ll>) J yO J 

( COSIpiuj) —SUllpiuj) \ 

I smil){uj) cosip{uj) j 

_ ( COS^ lp{u;) -sm2ip{u;)/2 \ 
\ -sin 211) (uj)/ 2 sin2^(w) )' 

At u = 0, the matrix elements of the octupole operator Qs have the following structure 



(A.3) 



where g is a finite real number. At 7^ the expectation value of is 

1(4)1 = |tr(qpH)| = \qsm2^ic.)\, (A.5) 

due to the assumption ii). 

At = matrix elements of the angular momentum operator are defined as 



ju,y jut, ] = f ^' 
3 ^lu jfifj. J V J 



(A.6) 

At 7^ 0, using the unitary transformation in Eq. (|A.1| ), we obtain for the matrix elements of j 



( ^■^^j'^; ^■^-j'^; ] = (A.7) 

/ jcos^V^M+j'sin^V^M (j'-j)sin2^H/2 \ 
1^ if - j) sin 2^{u;)/2 f cos2 ^(c.) + J sin^ ^(c) J ' 

From Eqs. (^) and it follows that KQg)] and \ji,ui^)\^ = (j' - j)^ sin^ 2V^(a;)/4 have a 

maximum at \ip{oj)\ = tt/A, where there is a strong mixing of the unperturbed eigenf unctions. 
Therefore, the fiu component of the Inglis-Beliaev moment of inertia 

2|Jm^HP ^ (/-jf sin^2V>M 

can be large enough to affect the total value of the moment of inertia due to the smallness of 
its denominator e^(w) + e^{uj) at this rotational frequency. 

sign of is not important for the absolute value of the octupole deformation. 
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Figure Captions 



Figure 1: The total binding energy in ^"^S as a function of the angular momentum. The 
calculated points indicated by symbols '+' and 'x' for ND and SD states, respectively, are 
connected by solid lines. Experimental values for states with different parities and spins are 
indicated by symbols: filled square is used for the positive parity and even values of /; open 
square is used for the negative parity and even values of /; filled triangle is used for the positive 
parity and odd values of J; inverse open triangle is used for the negative parity and odd values 



of /. The experimental values of the total binding energy and excited levels from Refs. ^5 
and IpBi) respectively. 



Figure 2: The angular momentum I versus the rotational frequency uj in ^^S" and ^^Ni. The 
calculated values of the angular momentum connected by solid lines are indicated by symbols: 
filled square and filled circle are used for ND and SD states, respectively, for ^^S; open square 
and open circle are used for ND and SD states, respectively, for ^^Ni; '+' and 'x' are used for 
the results of calculations with the Skyrme III interaction |]16[ . 



Figure 3: /^-deformations in ^^S* and ^^Ni as a function of the rotational frequency oo. The 
symbols are used here are defined in Fig.^ 
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Figure 4: The non-axial octupole deformation: (a) in the ND yrast band of ^"^S and (b) in the 
ND excited band of ^^Ni, - as a function of the rotational frequency lu. The calculated values 
connected by solid lines are indicated by symbols: filled square is used for |/53i|; open square is 
used for l/^ssj. SIII means the results of Ref. . 



Figure 5: Neutron single-particle energies near the Fermi surface versus the rotational frequency 
uj : (a) in ^"^S and (b) in ^^Ni. The occupied states with different simplex quantum number 
s connected by lines are indicated by symbols: filled square and open square are used for the 
s = +1 and s = — 1 states, respectively. The unoccupied states with the simplex s = +1 and 
s = — 1 are connected by solid and dotted lines, respectively. 



Figure 6: The density distribution: (a) in the ND yrast band of ^^S at u; = 2.2 [MeV/Zi]; (b) 
in the ND excited band of ^^Ni at u = 1.4 [MeV/;i]. 



Figure 7: The static, J'^^\lu), the dynamical, J^'^\uj), and the Inghs-Beliaev, Jmi^jj), moments 
of inertia versus the rotational frequency uj : (a) in the ND yrast band in ^^S"; (b) in the ND 
excited band in ^^Ni. See the definition of different moments of inertia in the text. 



Figure 8: Components of the Inglis-Beliaev, j7iB(w), moment of inertia as a function of the 
rotational frequency u : (a) in ^"^S and (b) in ^^Ni. Notations of components are defined in 
Eq. dg). 



Figure 9: The degree of mixing of intruder orbits, cos^-?/'(w), in the region of quasi-crossing: (a) 
in ^'^S and (b) in ^^Ni. The quantity cos^ il){uj) is defined in the text. 



Figure 10: The inverse of fluctuation of the octupole and quadrupole moments as a function of 
the rotational frequency uj : (a) in the ND yrast band in ^^5; (b) in the ND excited band in 



56 



Ni. The fluctuation AQi^ is defined by Eq.(|T^ 
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Figure 11: The potential energy curve along the jS^i direction: (a) in ^"^S and (b) in ^^Ni. The 
calculated values of the P31 indicated by symbols are connected by solid line. 



Figure 12: (a) The total binding energy and (b) the charge radius in nickel isotopes. 
Experimental values of total binding energies from Ref. A charge radius is calculated 

with the approximation = + .64 [fm^] from the calculated values of proton radius r^. 
Experimental values of Vc for ^^'^^Ni from Ref. for ^^Ni it is the extrapolated value from 
Ref. m. 



Figure 13: The total binding energy in ^^Ni as a function of the angular momentum. The 
calculated points indicated by symbols '+' and 'x' for SD and ND states, respectively, are 
connected by solid lines. Experimental values for different states are indicated by symbols: 
filled square is used for yrast non-rotating states; open square is used for the excited states 
whose binding energies are proportional to /(/ + 1) [h^] starting from the 0^ state; filled circle is 
used for other states. The calculations of the ground spherical state reproduce the experimental 
ground binding energy quite well. Experimental values of the total binding energy and excited 
levels from Refs. [p3] and respectively. 



Figure 14: Neutron single-particle energies of ^^Ni along the SD band near the Fermi surface. 
The symbols are used here are defined in Fig.^. 



Figure 15: The density distribution of SD minimum at = [MeV//i] in ^^Ni. 
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